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Abstract 

Using the LAMP model for nuclear quark structure, we calculate the binding energy and quark 



5^ structure of a i? meson with a D meson. The larger-than-nucleon masses of the two heavy quarks 

and the complete absence of quark Pauli Exclusion Principle repulsive effects produce a bound 
state with structure ranging from a deuteron-like state to a fully mixed four-quark state. Unlike 
the deuteron, however, pion exchange does not provide the dominant contribution to binding. 
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I. INTRODUCTION 

The relativistic Los Alamos Model Potential [H E] (LAMP) has been used to describe 
the binding and structure of ^He and ^He, including a good description [3] of the deep 
inelastic structure function of '^He. It does not, however, describe the deuteron due to the 
large separation of the nucleons and the dominance of single-pion exchange contributions 
there. The LAMP, lacking quark-exchange correlations, best encompasses medium and 
short-range meson exchanges (two-pion, p, and the like). It must therefore be supplemented 
with long-range single-pion-exchange contributions |1] for a better description of nuclear 
binding energies. 

In this paper, however, because the isodoublet B and D mesons are much more mas- 
sive than nucleons, the localization energy is much reduced. This brings them into closer 
proximity than the nucleons in a deuteron, or indeed, even in a large nucleus. Furthermore, 
the quark content does not require any pairs of quarks with the same (internal) quantum 
numbers, obviating the contribution of any of the (quark) Pauli exclusion effects that con- 
tribute to the short range repulsion between nucleons. Thus, this is a system in which one 
can expect greater accuracy of the LAMP and a significantly more deeply bound state than 
the deuteron. 

When this B-D bound state is observed, the deviation from our predictions here will 
provide a very good measure of the center of mass motion and breathing mode collective 
excitations. These are difficult to remove in the LAMP due to its relativistic nature. Since 
the non-relativistic analog model, the Quark Delocalization and Color Screening Model of 
Wang et al. [5] gives very similar results to the LAMP after removing such effects, we 
expect them to be small and hence, that our predictions here are reasonably accurate. 

A. Initial Concepts 

The LAMP treats the confining potential for quarks (and antiquarks) as a fixed scalar 
interaction in a Born-Oppenheimer-like picture, with the location of the potential minimum 
defining the system location. Quarks bound in a baryon or meson are treated as being 
bound within this potential rather than directly to each other. As such, there are immediate 
concerns about removing center-of-mass and breathing mode contributions to the evaluated 



state energy. This concern is ameliorated by comparing the energy of the interacting system 
of the two heavy mesons with the value at large (essentially infinite) separation. 

In this paper, in addition to the confining Lorentz scalar potential of the LAMP, we 
have included a Lorentz vector potential, as is required from the observed small spin-orbit 
interaction in the non-relativistic quark model. [B] In fact, the vector potential is also taken 
as linear, attractive, but without a Coulomb-like contribution, as discussed in Ref. [7]. 

In the LAMP, the confining potentials for each hadron are distributed in an array and 
are truncated on the mid-planes between them. While complex in general, for the case of 
interest here - two heavy mesons - the structure is very similar to that of the hydrogen 
molecule in the Born-Oppenheimer approximation, except for the linear vs. inverse distance 
form of the potential. In this case, the large masses of the c and b quarks further enhance the 
credibility of the approximation - they may be taken in the conventional heavy quark limit 
[S] as the fixed origins of the confining potentials for the light anti-quarks that complete 
each meson. 

At large separation between the heavy quarks, confinement guarantees the isolation of the 
light quark wave functions from each other. However, as the two mesons approach within a 
distance less than a few times their root-mean-square radii, the truncation of the confining 
potential allows for tunneling of each light anti-quark wave function into the confinement 
region of the other heavy quark rather than the one to which it is bound. The concept 
behind this is that a quark can only be confined to nearest center of color attraction, as 
in string-flip models [U], for example. This spreading out, or delocalization, of the wave 
functions naturally reduces the localization energy and provides an initial source of binding 
between the two hadrons. 

B. Color magnetic and quantum number issues 

In nuclei and other systems, this basic consideration is complicated by additional ele- 
ments: there are color-6 combinations of quarks and color-magnetic spin interactions of 
significance on the scale of the binding energy. Here again, the concerns raised by these con- 
siderations are considerably reduced - the color magnetic interactions between the heavy 
quarks are reduced to negligible values by their large masses. The light quark color mag- 
netic interactions with the heavy quarks are also significantly reduced. Only the light-quark 



to light-quark color magnetic interaction remains. This is at most ~ 50 MeV as seen [TO] 
in individual hadrons (nucleons, A's, light spin-0, and spin-1 mesons) where it depends on 
the color and spin-strong-isospin combinations determined by the constraints of statistics. 
Furthermore, here the presence of both color 6 and color 3 combinations, as well as spin-1 
and spin-0 elements, make it clear that strong cancellations of these color magnetic effects 
to low levels are to be expected. Therefore, in this paper we will only estimate these (in Sec. 
XI), ignoring details. We also will neglect the very small electro-magnetic contributions. 

Because of these simplifications, we will mostly ignore the fact that there are two neutral 
states {B^D^ and B^D^) that should exist and mix, splitting to form states of definite strong 
isospin (0 and 1) although both have J3 = 0. They also allow us to ignore the detailed spin 
structures, ranging from J = to J = 2, the latter with all of the quark spins aligned. Also 
unlike the individual nucleon case, the c and h quarks may combine anti-symmetrically to 
form a color 3 state or symmetrically to form a color 6. In the first case, the light anti-quarks 
must form a color 3 antisymmetrically, thus requiring the spin-isospin combination to be 
symmetric (1=0, J=0, or 1=1, J=l) and in the latter, the opposite is true - a color 6 and 
(1=1, j=o, or 1=0, J=l). 

Again, these allowed spin-isospin combinations for the light quarks would only produce 
significant energy differences if the color magnetic interaction were larger. The color 6 
combination of the heavy quarks would not be expected to produce any attraction, as 
indeed no such components appear in baryons, but the color 3 combination would. Neither 
of these effects is included here as the channel to color neutralization by decomposition 
into two color-singlet mesons {B and D) is almost open. In any event, symmetrization and 
antisymmetrization between the c and h quarks is moot as they are eminently quantum 
mechanically distinguishable. 

We turn now to the detailed calculations of the light (anti)quark wave functions in the 
double well defined by the Born- Oppenheimer- fixed heavy quarks. 

II. THE TWO- WELL WAVE FUNCTION 

For two wells separated by ± 5 at dimensionless [H] positions w± = {0, 0, ±5} along 
the 2;-axis (see Fig. 1), we define the wave function 

'^l{y) = ifjiyj) + eipiy^)^ where r-t = r + w-t = {x, j/, 2; ± 5} . (1) 



This represents, for example, a light M-quark (i.e., which we assume to be massless) mostly 
moving and confined in the well at r_ (the "left") provided by a heavy 6-quark. There 
may be some "leakage," represented by e, into the "right" well at r+, provided by, say, a 
heavy c-quark. We assume that the b and c quark masses are large enough to justify a 
Born-Oppenheimer approximation of this sort. There is a similar wave function "^li with r_ 
and r+ interchanged in Eq. (ll]) for a light J-quark mostly confined to the well at r+ with 
e- leakage into the well at r_. We wish to determine what the best values of the parameters 
6 and e are that provide a four-quark or molecular-like binding that form a bu cd system. 
The b and c are well separated compared with their Compton sizes. Since they have little, if 
any, wave function overlap and have distinct quantum numbers, anti-symmetrization issues 
are irrelevant. For the rest of the paper we will drop the subscript L on \1/, but it should be 
borne in mind when we finally compose the bu cd four-quark state. 

We take the t/^'s in Eq. ([I]) to be dimensionless four-component Dirac wave functions for 
light (here, assumed massless) u- and (i-quarks. They are solutions of 

Hd^P = [-ioi ■ V + V{t) + pSir)] ^j = E tjj . (2) 

Here V{r) is the time component of a Lorentz four-vector and S{r) is a Lorentz scalar 
potential (both to be specified below). With the Pauli spinor x assumed to be quantized 
along the z-direction with spin-projection rris, the normalized four-component s-wave Dirac 
wave function ^(r) is 

The upper and lower radial wave functions ipaif^) and ipb{f) can be chosen real. We have 
calculated them by solving the coupled radial Dirac equations [12] for (dimensionless) linear 
Lorentz vector and scalar potentials of the form [TT] 

V{r) = r — R and S{r) = r . (4) 

Here — -R is a negative displacement pushing the vector potential V{r) down below the scalar 
potential S{r), so that confinement trumps Klein-Gordon pair creation. [6] 

The solid curves in Fig. 2 show the calculated (dimensionless [H]) 15 wave functions ipa{r) 
and r%pb{r) when the potentials have R = 1.92, k^ = 1.253 GeV/fm. Physical dimensions 
are obtained by dividing the dimensionless r, R, etc., hj n = 2.52 fm~^ The ground 



state eigenenergy resulting from this calculation is 0.375 GeV. These potentials provide a 
reasonable fit to the c c spectrum. [7j 

For the calculations presented below, the ipa,b have both been fitted to sums of Gaussians, 

Mr) = E «^ e-'^''■'/^ Mr) = T.b^ e"^^^'/' , (5) 

i i 

where the Oj, bi, and /ij are dimensionless numbers. We found it necessary to go to six 
terms, so that evaluating the upper and lower components of the left-hand-side of Eq. ^ 
gives reasonable agreement with the right-hand-side. The fitted parameters are 

/ii= 1.0,1.3,1.6,2.0,4.0,8.0} 

ai = { 0.492649, -0.687482, 1.84609, -0.00246039, 0.258295, 0.0956581 } (6) 

bi = { -0.0571296,1.03367,-1.18398,1.33989,0.162575,0.299479 } . 

The fitted ipa{r) and ripb{r) are shown as the dashed curves in Fig. 2, largely overlying the 
solid curves from the solution of the Dirac equation. To check the quality of the fits we 
have evaluated the single-quark expectation value of the Hamiltonian, < H^, >= 0.7545, 
which is slightly larger than the (dimensionless) energy eigenvalue E = 0.7540 (which, for a 
variational trial function, is as it should be). 

III. EXPANDING (Hj^) 

The idea is that we will want to minimize the expectation value {H]j) (and the color- 
magnetic interaction between the quarks) with respect to the parameters e and 6 to bound 
(approximately) the energy for the four-quark system consisting of b, c, u, and d. {H^) is 
required for a variational bound as, due to negative energy states, < Hr, > itself is unbounded 
below. The Dirac Hamiltonian Hd is displayed in Eq. ([2]) but now, for the two-well case 
(Fig. 1), the potentials are 

r_ - i?, if z > _ r_, ifz>0 

V{r) = <^ and 5(r) = <^ . (7) 

I r+ - i?, if z < I r+, if z < 

As already mentioned, —R is a negative offset so the vector potential lies below the scalar. 
The exact two-weW energy E is in principle found by solving for the eigenvalue of 

Hj, ^(r) = E ^(r) , (8) 



with \1/ given in Eq. ([I]). This being difficult, we instead chose to find an approximate value 
of the four-quark energy E by the above-mentioned minimization of {H]j) 
After some algebra one finds 

Hi = -V^ + V\r) + S^r) + 2/3 V(r) 5(r) 

-ia ■ [{VV{r)) + /3 (VS(r))] - 2i V{r) aV . (9) 

The lack of a term like —2i S{r) ck ■ V is because of a cancellation as the Dirac operators 
CK and P anti-commute. The first four terms of H^ are "diagonal" (generically, Od) in 
that they connect tpa to ipa and tpb to ipb, while the last two terms are "off-diagonal" {Ood) 
connecting ipb to ipa- 

As a second check on our Gaussian fits of ipa and ipb, Eqs. ([s]) and (6), we evaluated the 
expectation < iJ|, > to be 0.5691, again slightly larger than E^ = 0.5685, as it should be. 

The following six sections describe how we calculate the expectation values of the terms in 
Eq. ^. These are basically integrals over various Gaussians, and are often messy formulas 
involving error functions. A reader less interested in these details may wish to skip directly 
to Sec. X, which presents the numerical results of these calculations. 

IV. GENERAL REMARKS ON CALCULATING THE EXPECTATIONS 

The reason for approximating our numerical radial wave functions ipa and ipb as sums of 
Gaussians is that it allows us to calculate the expectation values of each of the terms of H^ 
analytically. Given an analytic expression for H^ allows us to plot it quickly and precisely 
as a function of 6 and e. To do these integrations, we have relied heavily on programs such 
as Mathematica and Maple. As will be seen, the final results can sometimes be messy and 
often involve error functions [13j because of the Gaussians being integrated. 

For the diagonal operators of if|, we will calculate the upper and lower contributions 
separately, 

{^lOom = (^iOi,i^)^ + (^iodI^)s . (10) 

The 5-expectations are more complicated than those for A because of the factors of —icr ■ r± 
multiplying the radial ipbS. However, for some diagonal operators, as will be seen below, the 
5-expectations are not always needed. In any case, from Eq. ^ we expand these diagonal 



operator expectations as 

where the lij and Jij are integrals over Gaussians. 

First, we separate out the quadratic dependence on e as 






^I^^^I^ 



:l+e^)C+^4/^ 



(11) 



(12) 



and hkewise for the lower- component i?-integrals Jij. The second equality here conies about 



because parity symmetry ensures that the l\j = I^j , etc. We will refer to the I^j 



terms," in that they connect Gaussians with fij r^ /2 to those with fiir'^/2 (and similarly for 
Jj • with r+). Recalling the l/An from the normalization of the ip^s, we ensure the symmetry 
under the interchange of indices i and j by writing 

■^ 8n J ^ 



-IJ.irl/2 



(13) 



The direct integrals Jij have a similar form but with cr ■ r Od cr ■ r_ in place of the Od- 
The "cross terms" lij are more complicated integrals than the ^j , and likewise for J,-j- . 
They connect Gaussians with /ijr^/2 to fijr'^/2 and vice versa. Thus, on symmetrizing in 
i and j, 






- fd\ I 

Svr i ^ 



+ 






(14) 



The Jij have a similar form but with Od replaced by cr ■ r_ Od cr r^ or cr ■ r+ O^ cr ■ r 
as appropriate. 

Each of the off-diagonal operators in Eq. ^ has the general form 



O, 



OD 



-ioL ■ X 





-i cr ■ X 



21 



-i CT ■ Xi2 





(15) 



where the X12 and X21 are vector-operators that may not be equal because of the possible 
presence of the diagonal /3 matrix in Oqd- 

The direct terms of the off-diagonal expectation < Oqd > involve several terms because 
the upper component of \l/"'"(r_) connects through —i cr ■ X12 to the lower component of 
\l'(r_) at the same time that the lower component of \l/^(r_) connects through —i cr ■ X21 to 
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the upper component of \l/(r_). We therefore have to keep the sums over the a's and 6's in 
Eq. (Is]) as parts of the integrand. Again symmetrizing in i and j, 






(16) 



The cross terms of < Oqd > have even more terms because the \l/^(r+) connects to \E'(r_) 
at the same time that \E'"''(r_) connects to \l/(r_|_). It becomes 











+ ajbi{i 


+ e-/^.r2_/2 


-aibj{ 




+ ajbi(( 


+ e-/^^-+/2 


-ttjbii 




+ aibj{( 


^^-Hrl/2 


-ajbi{ 




+ aibj{ 



cr-Xi2)(cr-r+) 
a-r_)(o--X2i)] e-^^-^+/2 
cr-Xi2)(cr-r_) 
^•r+)(^.X2i)] e-^»^^/2 
(T ■ Xi2)(cr ■ r+) 
o--r_)(o-.X2i)] e-^^'-+/2} 



(17) 



The integrations for the J's, J's, and in Eqs. (16) and (17) can best be done using (dimen- 

1 /9 

sionless) cyhndrical coordinates, p = (x^ + y'^) , 6, and z. The 6 integrations are trivial, 
providing a factor of 2tt, which will cancel with the I/Att coming from the normalizations 
of the ip^s in Eq. ^ to give an overall factor of 1/2 before each double integral over p and 
z. It usually is easier to do the p-integration (from to +oo) first. Because V"(r) and 5'(r) 
depend on r_ when z > and on r+ when 2; < 0, we need to do the ^-integration separately 
for those regions, i.e., for z from —00 to and then for z from to +00. The separate 
results are then added and simplified to give the final integral. 

We will distinguish the results for the expectations of the different operators in Eq. ^ 
by an appropriate subscript. For example, for Od = V^, we will write Jj-^' as -^i'<;v2>' ^^^ 
similarly for the Jij integrals. 



V. NORMALIZING ^ 

While the Dirac ip^s are themselves properly normalized, the two- well "^ is not. For this 
we need to calculate the expectation values of Od = 1 to find 



N^{5,e) = j d^r^^m = (^|1|^) = (^|1|^)^ + (^|1|^) 



B 



(18) 



We make the expansion in e as in Eq. ( 12 ) above. The direct-term integrals for the expec- 



r(0) 



tation ( 1 ) are, noting that for the Jij'^iy we also have a factor of (cr ■ r_)(cr ■ r_) = r^ in 
the integrand. 






71 



nl/2 



_2{fli +/ij)^_ 



t(o) _ 3 
-'jj, <1> - -3 



TC 



2(/ii + jjj) 



A5 



1/2 



both independent of 6. 



The cross-term integrals do depend on 6. For the ^jj <i> we need the factor 



(cr • r+)(cr • r_) = r+ ■ r_ = p + z —6 



(19) 
(20) 



(21) 



in the integrand. Proceeding as in Sec. IV, we find 



ij, <1> 



27r 



.1/2 



-2ij.iHj S'^/ifii+fi.j) 



J, 



(1) 



ij, <1> 



_(/ii + /ij)3_ 

3(/ii + /ij) - 4 /ij/ij (5^ 



27r 



(/^i + /^i)' 



1/2 



'2fiifij 5'^ /{^li+^ij) 



r(l) 



(22) 
(23) 



Note that, when 5 = 0, -^ij/<i> = 2 I^j'^i^, and Jjj/<i> = 2 Jjj;<i>. This is a common 
feature for all the expectations here and below. This is necessary so that, for example, when 
6 = and e = 1, one recovers a result that is four times that when 6 = and e = 0. 



r(l) 



We see from Eq. (22) that Jj,- <i>,as a function of 6, is a decaying Gaussian (as in Fig. 



SK). On the other hand, Jjj;<i> falls off from its peak at 6 = 0, goes through zero, and has 
a mild minimum before decaying to zero at large 6'^ (as in Fig. |3| 
Combining all terms. 



N'{e,6) = j:a.a, (1 + 6^) 41 + e4V> + E^^ &. (l + ^^^iV + ^^V 



(1) 



.2A 7(0) 



r(i) 



(24) 



«j 



«j 
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and the normalized \1/ is obtained by multiplying Eq. ([I]) by 1/A^(e, S). 

Figure [4] shows a plot of A^^(e, 5) for the values of the a's, 6's, and /i's that were fitted to the 
normalized ipa and ^pb, Eq. (|6|. We have checked that, for these values, A^^(0, 0) = 0.9858 ~ 1 
and A^^(1,0) = 3.9430 ~ 4, as they should but with some deviation {^ 2%) coming from 
the inexactness of the fitting. The ratio of the two values is 4 to high accuracy. 

To illustrate what "leakage" from one well to the other might look like. Fig. [5] shows a 
plot of the upper component of the normalized \1/ as a function of p (running from to 2) 
and z (running from -3.5 to +3.5) for e = 0.5 and 6 = 1.1. 

VI. EVALUATING THE DIAGONAL EXPECTATION ( - V^ ) 



First, note that, for r± = {x, y, z ± 6}, the ith component of the gradient 

d , d 

Vj = — = V- = — for r' = {x = x, y' = y, z' = z ± 6} = r± 



(25) 



since each dx'Jdxi = 1. Thus we can replace the result of the Laplacian with respect to 
r acting on a function such as ipai^-) with that for a Laplacian with respect to r_ acting 
on that function. For spherical coordinates, — V^ on the angle-independent e"'^^'"-'^ then 
becomes 

1 £ 



-V^ e"''^' -/^ = -V'^ e"''^ -/^ 



r_ dr"^ 



( r_e-'^^'^-/2 ) = -^^^^.rl-3) g-^^''-/^ , (26) 



whence the three-dimensional integral reduces, after symmetrizing and canceling factors of 

1/2 



47r, to 



7-(0) _ q 

-'ij, <-v2> ~ '^ ^^i^^j 



IT 



2{ni + Hj)^ 



(27) 



independent of 6. 

For the i?-integrals things are more complicated because of the cr ■ r_ factor to the right 
of the Laplacian. After some algebra. 



(cr ■ r_ VV ■ r_) e-^"^''-/^ = -rl fi^ {p^rl - 5) e-'^^''-/^ 



(28) 



whence 



also independent of 5. 



t(0) 



15 ^ii^ij 



TT 



2{fii + pj] 



1/2 



(29) 
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The cross terms, again, do depend on S. 



y-(l) 



27r 



.1/2 



_(/ij + ^ijY _ 



-2ni/ijS^/{fii+fij) 



(30) 



A plot of this integral as a function of 6 looks like Fig. |3j 

For the corresponding i?-cross term, one proceeds in the same manner but, instead of 



Eq. (28), we need [H] 



{a- ■ r+VV ■ r_) e''^^'--/^ = -(p^ + ^2 _ ^2) ^^ ^^^^2_ _ 5) g-M,.2_/2 



(31) 



We find 



r(i) 



2, .2 r4- 



•^i, <-v2> = AiiAij [ 15 (/Ui + /Uj) - 40 fiifijifii + fij) 5+16 fi^fij 6 



X 



27r 



(/ii + /ij)ii 



1/2 



-2fj,iij.j5^/{fii+fij) 



(32) 



This integral as a function of 5 also looks like Fig. |3]B, but because it is quartic, it is slightly 
positive beyond 6 = 1.7. 



VII. EVALUATING THE EXPECTATION OF F^ + 52 _^ 2/3 VS 

This is also a diagonal operator. The linear vector potential V"(r) differs from the linear 
scalar potential ^(r) by a negative offset — R. In (V"^(r)) the integrals of (^±) are the same 
as those for (^^(r)). Here, (^±) means the integration of r^ when z > and of r^ when 
z < 0. Thus we (schematically) expand the diagonal V"^ + S^ + 2(5 VS as 



V' + S' + 2(3VS) = 2(ri){l + l3)-2R{r±){l + l3)+R'{l). (33) 



The factor of (1 + /3) ensures that only the upper components of ^ contribute to the first 
two expectation values. That is, we only need to calculate the A-integrals (the J's) for 
those terms. The expectation value ( 1 ) multiplying R^ does have contributions from the 



lower components and their integrals /<i>, -?^<i>, J<iyi and J<i> are given in Sec. 



V 



The 



integrals for the operators (^±) and {r±) are rather more complicated and their analytic 
forms are presented next. 
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A. Od = t\ 



The direct integral for this operator is 



7-(0) 



26 



if^i + f^jf 



--{^l^+^lJ) 52/2 



+ 



TT 


1/2 


2{^^i + ^^j)\ 



3 + 2(/x. + /.,)<52Erfch/^^'^^5 



(34) 



Note the hnear dependence on 5, which gives rise to a shallow minimum near the origin 
before the function returns to its initial value, as in Fig. [3p. 

The cross-term integral for < r^. > is 



J ■ ■ ^ 2 ^ 

^3, <r-^> 



4(5 



+ 



27r 



-(^ii+Atj) 52/2 



1/2 



-2^li^JLj 5^/{tii+tij) 



3(/i, + /i,) + 2 (/i^ + ^2) 52 _ 2 (^2 _ ^2) ^2 E^f 



(35) 






which also has odd terms in S. In this case, as a function of S, l\- ^^.2 > falls off smoothly to 
zero from its peak value at 5 = 0, as in Fig. 30. /• ,^2 > is symmetric in i and j because 
Erf(-x) = -Erf(x), /^^^<^,|> = ^]J^<r|>- ^^^o, as expected. 



j(i)^ = 2 /(°l = 3 



2tx 



(/ii + /ij)^_ 



1/2 



(36) 



when (5 = 0. 
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B. Od = r± 



The direct term for this operator is 



r(0) 



1 



rV"; _ ' [4-1-9 p-2(M»+A'j) <5^ _ Q p~{t^r+H) <5V2l 



1 

25 



TT 



2(/ii + /ij)5 



1/2 



(37) 



(/ii + /ij) 5^ - (l + 4(/ii + /ij) (5^) Erf (^y^2 (/x, + ^j) 5^ 
+ (l + 3(^, + ^,)52) Erf(^v'(/i, + /i,)/2 5)| . 



This integral also has an odd term in 5, like I^- 2~^- Its plot as a function of S resembles 
that shown in Fig. SC. That is, despite the 1/5 factor in the last term, l\j^^^^y is not singular 
at (5 = (i.e., when there is no separation between the two wells): l\j.^j.^^ — )■ 2/(yUj + /ij)^ 
as (5 — )■ . 

The cross term for {r±) is 



'■I], <r±> 



1 

2 (5 /ij/ij 



/'g-2M, 52 _^ ^_2^^ 52 _ ^_l(^^+^^,) 52N, 



TT 



2(/ii + /ij)^ 



1/2 



(/.. + /.,)^ErfcfJ^^5 



2i5 /j^pj 



- 2 /ij (^/ii + fij + A fi^ 5^) e ''«+^j Erfc 



MjMj 



+ 2 /ii (/ii + /ij + 4 /i| 5^) e ^^+^^ Erf 






/ii5 



yUi + yUj 



/ij5 



2<5 fjjfjj- 



(/ij — /ij) f/ij + /ij — 4 /ij/ij (5^ j e '"'"''''j Erfc 



^^2 {fii + fij) ^ 



(3^ 



Note that L 



(1) 



we have /, 



(1) 



is also symmetric under the interchange of i and j and, again, at 6 = 0, 
- 4:/{fii + fij)'^ = 2 I^j^^j.^^. Its behavior as a function of 5 is similar to 



«J: <r±> 

that shown in Fig. [3p, again partly due to the presence of odd terms in 6 



VIII. THE OFF-DIAGONAL EXPECTATION OF -ia [(VF(r)) + /3 (V5(r))] 



For the linear potentials of Eq. ([T]) 

V\/(r) = V5(r) 



f_ if 2; > 

f + if 2; < 



(39) 
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and we again have a simplification from tlie (1 + /3), namely, 



za ■ [(VV(r)) + /3 (V5(r))] = -ta ■ f±(l + /3) 





-2i cr • f± 



(40) 



i.e., the operator X12 in Eq. (15) vanishes and X21 is doubled. The latter operator connects 



the upper component of \1/^ to the lower component of \1/. 



For the direct terms, Eq. (16) reduces to two terms 



< [VVSf^ > = -'^^Y.Jd'r [e~^^'-/^ [a,bi {a ■ r_)(cr ■ f±)] e"'^^"-/' 



where 



^gUy5> = -2 ^ Jd'r e-'^-^/^ [{o- . r_)(^ ■ f±)] e" 



^ij,<S7VS> 

The Pauli matrices here reduce to 



-Mj'--/2 



[cr • r_j (cr ■ r±j = r_ (r_ ■ r±, 



(41) 



(42) 



(43) 



For the integration over z > the integrand becomes simply r_, which is the same as that 

above. For the 



(0) 



already needed for getting to the final result for /<r.j^> in subsection 



VII B 



integration over negative z, however, Eq. (43) becomes 



r_ r_ ■ r I 



r_ ■ r+/r+ = {p' + z' - 6')/Jp^ + {z + 6) 



(44) 



which involves a new integrand, but which nonetheless can still be done analytically. (Here 
it is much easier to do the p-integration first.) We find 

2 



K. 



(0) 
ij,<\7VS> 



{pi + Pj] 



1 

~5 



27r 



2 — g-C^i+Mj) -5^2 
1/2 



(Pi + /ij)5 



Erf ^p{pi + /i,) 5 j - Erf (^^(/i, + p,)/2 5 



(0) 



(45) 



This result is, again, symmetric and non-singular with -^ij/<vy5> = ~4/(/ij + Pj) at 5 = 0. 
In this case there are no odd terms (!) in 5. Its plot versus 5 is similar to that shown in Fig. 
[3p, but with the initial slope at the origin being zero. 
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Because Kij^<^^vs> — ^ji,<vvs>i ^^ ^^^ finally write the direct term contributions for 
this expectation as 



< [VVSf'^ > 



y . ('^i"« + ^i'^j) ^ij,<\/VS> ' 



(46) 



«j 



regaining explicit symmetry. 

The cross term integral if^yy^^ is more complicated but is done similarly. As X12 = 0, 



there are now four terms remaining from Eq. (17) 



< [VVS]^^'^ > = -2 — E /^'^ {e-^^-+/2 [a^hi (cr ■ r+)(cr ■ f±)] e 



IMjr'i/2 



+ e 



-Ai,r2/2 



K-6,(o-.r_)(o-.f±)] e-^^'+ 



H d/^ 



Y. d 



«j 



ttjOi -f^ij,<vy5> + '^i'^j ^ji,<vvs> 



(47) 



where 



K. 



(1) 

ij, <vvs> 



— i d'r it 
47r 7 1^ 



^tJi^rl/2 



[(o--r+)(o--f±)] e 



-/.,r2/2 



+ e 



-^^^r1|2 



[(o--r_)(o--f±)] e 



-A'r'-V2 



(48) 



In addition to Eq. (43) we also need 



{a ■ r+j (o- • r±j = r+ (r+ • r±, 



(49) 



which becomes r+ for the 2 < integration and (p^ + z^ — 5"^)! J p^ + {z — 5Y for the z > 
integration. 

The integrations over z go much easier if one re-defines the integrations over z in terms 
of /i = /ij + /ij and u = fii — fij. The resulting integrals in /i and z/ can then be converted 
back to fii and fij. We find 
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K 



(1) 

ij, < VV'5> 



fiifijifii + fijf 



2 /ij (/xj - fii) e-2^' ^' + 2 /ii (/ii - fij) e-^^^ ^' 



1 



2 5/if/i 



2, ,2 



TT 



2{fii + fij] 



1/2 
X 



2 ^-(m>+Mj) 52/2 



(50) 



|(/Xi + fijf (/Ui + /Uj - 2/ij/ij (5^ ) Erfc f ^ (/ij + /ij)/2 5 
+ 2 /i,^ {fi^i + Afiifij + 3/i^) - 4/i^(/ij - /ij) (5^ 

X e-^^^^"^ ^V(P^+H) Erf ('y'2/(/ii + fij) fi^ 6 

- 2 /i^ [ 3/i- + 4/ii/ij + /i| - 4/i,2(/ij - /ii) 5^ 

X e-2'^''^^ ^V(M.+/^.) Erfc ('y^2/(/i, + /i^) fn 5 

- (/^i + 5 /i • /ij + 5 /ii/i| + fi^j) - 8 /i^/i| 6^ 

X (/i, - /i,) e-^^'''^) '5V{M.+/^.) Erfc | i^fi^Ml 

VV2(/ii + /ij) 

which also is symmetric and goes to — 8/(/ii + /i^)^ = 2 K^j^l^^yg^ ai 6 = 0. This integral 
does have some odd terms in 6. Its plot as a function of 6 resembles a Gaussian, i.e., looks 
like that shown in Fig. [3]A.. 



Because -f^ij/<vys> — ^ji,<vvs> ^^ "^^^ again finally write 



(1) 



< [vvsr^ > = y: {a,h + a,b,) Ki;;^ys^ 



(51) 



ij 



mirroring the form of Eq. ( 46 ) . 



IX. THE OFF-DIAGONAL EXPECTATION -2i V{r) ex V 



For this off-diagonal operator X12 = X21 = —2V{v) V in Eq. (15) and the direct term 



expectation, Eq. (16), has all four terms 



«j 



<[2VVf^> = ^ Y,fd^rV{r) {g-'^'^'/^ [ 2a,6^.(cr ■ V)(cr ■ r_) 

-2biaj{(7-r^){(7-V)] e-^'^'-/^ 
^^-Hrl/2 [2aj6,(cr. V)(cr-r_) 

-2b,a,{(T • r_)(cr ■ V) ] e"'^''--/'} . (52) 
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With 



-/i, 



(r_)fce-^»^-/^ Vfc(r.), = 4;, and (o--V)(o--r_ 



(53) 



we have, for the first terms in the square brackets of Eq. (52), 

(cr ■ V) (cr ■ r_) g-'^'^-/^ = e"^'"-/' (cr ■ V)(cr ■ r_) + cr ■ [(cr • r_)V e''^'^-/^] 
= [3 - /x,(o- ■ r_)(a . r_)] e-'^-^-/^ = (3 - /i, r^) e''^^'--/^ 

and similarly when acting on e"^-'^-'^. 

For the second terms in the square brackets of Eq. (|52| , 



(54) 



(cr-r_)(cr- V) e 



-Miri/2 



-/ii(cr ■r_)(cr ■ r_) e^ 



-Mi'"i/2 



-/iir^e-'^^^-/^ 



(55) 



and, again, similarly when acting on e 



-fi,rl/2 



With Eqs. (54) and (55), Eq. (52) reduces to 



id 



-M,r2/2| 



E {«^&. < 



<2VV> ^^ "j"« -"-jj, <2yv>/ 



(56) 



«j 



where, with V"(r-|-) = r± — R, 



kS!<2VV> = 1, Jd'r e-^^ "-/^ (r± - i?) [(/.. - ^,) r^ - 3] e" 



.(0) 



'(0) 



.(0) 



Hrl/2 



'(0) 



= (^, - /i, ) i^; - (/i, - /.,) i? K^; + 3 i^:, - 3 i? i^:, 

where these four integrals are 

'(0) 



(57) 



^' Air 



2 (/ij + /ii 



+ 



26 



vr 



Ig _l_ Q g-2(/ij+/ii) (J^ _ ]_]_ g-(Mj+Mi) -J^/a 

.1/2 



2(/ij +/ij)'^_ 



[5 + 9{i2, + fii)S'] Erfc J(/i, + /ii)/2 5 



[5 + 12(/^, + /^i)52] Erfc (^v'2(/i, + /i,) (5^ | 



40) ^1_ / ^3^ g~^,ri/2 ^2 g-;.,.^_/2 ^ 3 

^' 47r 



7r 



4:71 



Ko) 



2(/ij +/ij)^ 



1/2 



Kgl, = ^ydVe---/^e-^^-/^ = /gU. 



(58) 
(59) 

(60) 

(61) 
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KIj[^ has an odd term in 5 and its plot resembles that shown in Fig. pp. All four of the 
above integrals are symmetric in i and j, so we can finally write 



(0) 



'(0) 



< [2VVf^ >=Y: ia,h + a,b,) K^;^vv> = 2 E«.&^ ^S;<2yv> • 



(62) 



«j 



For the cross term, from Eqs. (54) and (55) and the like, Eq. (17) becomes 



'^ hi 
{ e-^^^y^ [2a,hj{(T ■ V)(cr ■ r_) - 26^0,(0- ■ r+)(cr ■ V) ] e-^^'-'^ 

_^g-M.ri/2 [2aihj{(T-V){(T-r+)- 2hiaj{(T-r_){(T-V)] e-^^'H'' 

+ g-M.r-y2 [ 2a,-6i(cr ■ V)(cr ■ r_) - 26,-ai(cr ■ r+)(cr ■ V) ] e"^'"-/' 

^g-M,r-^/2 [2a^.6,(cr-V)(cr-r+)- 2 6,ai(cr ■ r_)(cr ■ V) ] e-^'^'+Z^ | 






-Mjr?./2 



+ e" 



+ e 



+ e" 



-^,r^/2 



-M,r2/2 



-M,r2/2 



aibj{3 — /Uj-r^) + hittjUjir^ • r_) 
aj6j(3 - /ijr^ ) + bjaifii{r+ ■ r_) 
aj6i(3 - /ijr^) + bjaifii{r+ ■ r_) 



-/.,r2/2 



-Miri/2 



7 . I '^i"i -^ii,<2VV> + S"» ^ji,<2V\/> / 



(63) 



«j 



where 



K, 



(1) 

jj, <2 y v> 



47r 



d^r (r± - i?) { e-'^'^+Z' [(3 - /i,r!) + /ii(r+ ■ r_ 



,-Mjr^/2 



+ e- 



-fj.irl/2 



(3- fijrl) + fXi{r+-r_] 



o-H^l/^ 



-fx,Kij^, + l^jRKlj\ + fi,Ki^,-fx,RK^, + 3K^Jl-3RKff^. 



(64) 



The first integral, 

K^i =— I dh { e-^' ^+/2 ^^ ^2 g-^, r2_/2 ^ g-^, r2_/2 ^^ 

"'' 47r J '- 



r^e-^^'-+/2} 



(65) 



can be done using /i = /ij + /ij and z/ = /i, — /ij, noting that ;U > |z/|. Writing 



g-M,r2/2^2 g-A«jri/2 _^g-Miri/2^2 g-M,-r2/2 



2 e-'^('''+^'+^')/2 { (p2 + ^2 ^ ^2^ cosh(z/fc) + {2z5) smh{v5z)] (66) 
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displays the i, j symmetric and anti-symmetric parts explicitly. After converting back to /ij 
and fij, 



K. 



(1) 

ij,a 



H'jKH'i ' H'j / 



~2iJ.i <52 



- { [ 5 /ij {fii + fij) +4:fi^fij 6^ ] e 



+ 



2 5 iXiiji] 



TT 



2{^Ji^ + ^ljY 



1/2 



(67) 



I 2 [ /ij (/ij + /ij)^(2/ij + 5/ij) + Ajj,ijj,j{jj,i + Hj){i4 - 2^i^j + 3/i^) 6^ + 16 /i^/i^ (5^ 

^ g-2M./., 5V(M.+M.) Erf ('y^2/(/ii + /i,) /i,- 5 \ 



4 x4 



- 2 ^;[ 3 (/ii + /i,)' + 24 iii{iii + /i,) 5' + 16 /i^ 5 



^ g-2M./., &V{^^^+N) Erfc ('y^2/(/ii + /i^) /i^ 5 



2\ x2 



- 16 /i^/i^ (/i, - ^^) 6'] e-^/^'^^ '5V{A«.+/^.) Erfc I ^^^ ~ ^^'^ 5 



+ {fXi + fXjf (2 /ii + 3 /ij) Erfc (^y^(/ii + /ij)/2 5 j | 



2(//i + /ij) 



which is, as expected, not symmetric in i and j. It is, however, non-singular: K, 



(1) 



16/(/ij + fij)^ at (5 = 0. Its plot resembles that in Fig. |3]D. 
The second integral is much simpler. 



^'^'' 47r 



hh'^{ 



27r 



g-M»'-+/2 ^2 g-Mjr-?,/2 _^g-M»r-2/2 ^2 g-Mjr-2_/2| 
1/2 



(yUj +/ij)'^ 



3 (/ij + /ij) + Afi^j 5 



-2fiifj,j5^/{Hi+^ij) 



(68) 



which is also non- symmetric, but only because of the term proportional to S"^. As a function 
of 6 it looks like Fig. |3l 
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Almost as complicated as KI,\^, the third integral is 



K, 



(1) 



Att J ^ 



+ 2/ii [ (/ii + /ij)(4/ij - ^i) - 4/ii/i2 52 ] e-^'"^ ^' 
+ [ (;Ui + /ij)(/i- - 8/ii/ij + /i^) + 8/i-/i^ (5^ ] e" 



;(m»+Mj)<5^ 



+ 



2 5/.f/i2 



TT 



2(/ii + /ij)' 



1/2 



(69) 



I 2 /i^[ (/ii + /ij)^(4/ii + /ij) + 8 /i,^ (yUj + yUj)(2/ii - /ij) 5^ - 16 yuf/ij 5"^ ] 

^ g-2M,A', 5V(/.,+/^,) Erf (J — ^ — /ii 5 1 
- 2/i^[ (/ij + lijfiyAiij + /ij) + 8/i^ (yUj + fij){2fij - yUj) (5^ - 16 /ii/i| 5^ 

^ g-2M,Mj 5V(/^.+/^.) Erfc I J — - — a. 5 ] 



x(/ii-/ij-)e'^'''^^^'/^^'^''^^ 



+ (/Xi + /ij)^ [ (/i- + 3 /ii/ij + fl'^j) - 2 yUi/ij (/ii + /ij) 5^ 



1 + Erf( 4^i^i^5 
2(/ii + /ij) 



X Erfc (^y/(^, + ^,)/2 5)1 



Ki) 



(1) 



'(1) 



which is surprisingly both symmetric, -fQj/c = -^ij/c ^^^ non-singular: K^j'^ = 16/(yU.j + /ij 



at (5 = 0. This integral as a function of 5 looks like Fig. |3p. 
The fourth integral is also simple, 






ih'^i 



27r 



1/2 



(r+-r_)e-^'''+/2| 



(/ij +/ii)^. 



[ 3(/ij + /ii) - 4/ii/i,- 52 ] e-^^'^^^^ -^ /(^'+'^^) . (70) 

Its 5 dependence. Fig. |3]F, shows a relatively deeper minimum than that depicted in Fig. 
3p. The fifth and sixth integrals are already familiar, 

(1) _ L f d^r (e-'^^^'V^ ^. .-M,r-i/2 , „-M»ri/2 ^. ^-^,r-i/2\ _ jW 

[it J ^ 



K. 



«J,e 



47r 



^^e-^^^-/2 + e-^»''-/2 r±e-^^^+/2} 



'ij, <r±> 



(71) 









-M,r2/2 „-M,r2/2 



J U 



(1) 



<1> 



(72) 
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These last three integrals, -ft'L rf through K^, K, are all symmetric in i and j. 



X. PLOTTING Hjj ^ TO FIND A MINIMUM ENERGY 

We can now pull all these integrals together to get an analytic expression for < Hj^^ >, 
from which we can look for a minimum squared energy. First, we define the (unnormalized) 
contribution, as a function of e and 6, from the diagonal pieces, 

Hn%ae, S) = Y.a. a, [ (1 + e') (4°^ ,, + 4 1^^^^^^^ - 4 R I^^,^^ + R' l\f,,^) 



«j 



+ e ( h],<V'> + 4 h],<rl> - 4 i? Iil<r±> + R^ hl<l> 



+ E &. &, [ (1 + e') (4°!cv^> + R' 4;^i>) + ^ (4;^v^> + ^^ 4;^i>) ] • (73) 



h3 



Figure p^ displays a three-dimensional plot of Hj^'^^^^J^e,5)/N'^{e,5). It shows a relatively 
shallow minimum at e = 1 and 5 ~ 0.8. Note the large value, a dimensionless squared- 
energy of ~ 4, which must be largely canceled by the off-diagonal contributions to achieve 
a squared-energy similar to that for the one- well case, E"^ = 0.5685. 
The off-diagonal (unnormalized) contributions are 



-"D, off-diag(^' o) — 2^ Oj bj (1 + e ) ( -f^jj,<vys> + ^ij,<V^> 



+ e (k^^<vvs> + <<yv>) 1 • (74) 



Figure M gives the three-dimensional plot of Hj^'^^^_^^^J^e,5)/N'^{e,5). In contrast with 
H o'^dia.g/ ^"^ 1 i^ has a repulsive hump around (5 ?« 1 as well as a shallow valley running 
from e = to 1 for 5 ~ 0.2. In the final sum of diagonal and off-diagonal contributions that 
hump will fill in the minimum seen in Fig. |6j 

Thus we finally combine the two contributions, defining a normalized 

H^\e, 6) = [i/^,%_diag(e, 5) + i/z5,'off_diag(e, 5)] /N\e, 5) . (75) 

Figure [s] plots how -ff^)^, as a function of e and 6, develops a long, flat valley for all values 
of e at a separation of 5 ~ 0.2 (i.e., recalling the value of k, a separation of ~ 0.45 fm). Also 
important is the hump (reminiscent of a fission barrier) around 6 ~ 0.9 that will help to 
confine this four-quark system at 5 ~ 0.2. This hump corresponds to a repulsion between 
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two Q — q asymptotic meson states if the light quarks delocahze. There is very httle, if any, 
barrier to coalescence at e = 0. 

It is easier to see this behavior with a two-dimensional plot, Fig. ^ showing Hj-,'^ as a 
function of 5 at e = 1, where the valley is deepest and the hump is highest. 

The dimensionless squared-energy valley-depth at 5 = 0.18, Aif^^ = 0.097, corresponds 
to a binding energy of 155 MeV for this hcud four-quark mesonic state. The valley is 



surprisingly flat, as shown in Fig. 10, dropping only 0.0023 from e = to e = 1, correspond- 
ing to an energy drop of about 24 MeV. This suggests that Zitterbewegung may play an 
important role in the nature of this meson. 

XI. DISCUSSION 



Figure [TT] is a contour plot of the binding energy of the state in the e-5 plane. It displays 
two remarkable features: The first is that, at very small e, appropriate to the approach 
towards each other of the two asymptotic {B and D) mesons, there is no evidence of a 
repulsive barrier to the fusion of those mesons. The second is that the valley of attraction 
at small meson separation is very fiat between small e (~ 0.2) and e = 1. This indicates 
that there is little energy associated with fiuctuations in the e collective variable of the 
light quarks in the state. There may be a more significant amount associated with the 5 
collective variable, but this effect is suppressed by the large masses associated with the Born- 
Oppenheimer centers defined by the heavy quarks, at least when viewed non-relativistically 
as seems appropriate for them, due to their relatively large masses. We therefore expect 
little correction to our estimates of the mass of the four-quark state due to collective variable 
effects. 



The dashed curve in Fig. 11 illustrates how two well-separated B and D mesons at e = 
and large 5 would come together to 5 ~ 0.2, corresponding to a heavy quark separation of 
K. 0.45 fm. This small separation makes it clear that long-range pion-exchange effects can 
be expected to be negligible. From e ~ 0.2, the four-quark state then slides down the nearly 
fiat valley to e = 1 where it is most bound. Such a state is prevented from falling apart 
because of the "fission barrier" around 5 ~ 0.9. 

We have ignored the possible color magnetic contributions from the interaction of the 
two light antiquarks, but this must be less than 50 MeV and we expect it to be even less 
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than half this value. This correction is not large compared to the extracted binding energy 
of order 150 MeV found in our calculations. Thus, by comparing our binding energy with 
the threshold for B and D mesons, we predict a (set of) states in the region of 7 GeV/c^. 

If our calculation may also be applied to a pair of c and c quarks, there would be a four- 
quark bound state around 3.6 GeV/c^. Although exotic, the recently discussed X(3872) with 
J^'-^ = l"*"*" [15] is too massive to correspond to any of the states implied here. However, if it 
included a strange quark instead one of the u or d quarks that we have explicitly considered, 
and we compare to the mass difference between the D and Ds mesons, it is conceivable 
that it would match with our calculations here. Moreover, as is well known from hadron 
spectroscopy and lattice studies, the light mass of the s quark does not seem to have a 
significantly different wave function or a different behavior from the lighter quarks, so a 
recalculation along the lines we have described here might well be relevant to this state. 
(Again, the issues of symmetrization vs. antisymmetrization would be completely obviated 
as all of the quarks would be distinguishable.) 

Finally, we comment on the surprisingly small difference in binding energy between the 
"molecular" form of the bound state, (e ~ 0.2, as in nuclei [1])) and the four-quark limit 
(e = 1). If this feature is widespread in such heavy quark systems, it could go far towards 
explaining why it has been so difficult to identify unambiguous four-quark states. 
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FIG. 1. Two- well linear potential. 




FIG. 2. Massless quark Is wave functions ipair) (above the axis) and ripb{r) (below). 
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FIG. 3. Typical plots of /'s, J's, and K^s as functions of 6. 




FIG. 4. Three-dimensional plot of N'^{e,6). 
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FIG. 5. Plot of ^aip, z) for e = 0.5 and 5 = 1.0. 




FIG. 6. Plot of the diagonal contribution to H^'^{e,5). 
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FIG. 7. Plot of the off-diagonal contribution to Hj^'^{e,d). 




FIG. 8. Plot of the final H^'^{e,6). 
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FIG. 9. Two-dimensional plot of H^^{e = 1,6), showing the valley at 5 = 0.18 and the "fission 
barrier" at 6 ~ 0.9. 
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FIG. 10. Plot of how the nearly flat valley at 6 = 0.18 decreases from e = to e = 1 
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FIG. 11. Contour plot of Hjy'^{e,6). The dashed curve illustrates how two well-separated Q — q 
mesons at e = and large 6 would come together and slide down the valley at (5 ~ 0.2 to form a 
four-quark state at e = 1. 
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